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SUMMARY 

The deployment of large shares of stochastic renewable energy, e.g., wind power, may bring impor¬ 
tant economic and environmental benefits to the power system. Nonetheless, their efficient integration 
depends on the ability of the power system to cope with their inherent variability and the uncertainty 
arising from their partial predictability. Considering that the existing setup of the European electricity 
markets promotes the spatial coordination of neighbouring power systems only on the day-ahead mar¬ 
ket stage, regional system operators have to rely mainly on their internal balancing resources in order 
to guarantee system security. However, as power systems are forced to operate closer to their techni¬ 
cal limits, where flexible generation resources become scarce, the conventional market paradigm may 
not be able to respond effectively on the wide range of uncertainty. 

The operational flexibility of the power system depends both on the technical parameters of its 
components, i.e., generators and transmission infrastructure, as well as on the operational practices 
that make optimal use of the available assets. This work focuses on alternative market designs that 
enable sharing of cross-border balancing resources between adjacent power systems through High 
Voltage Direct Current (HVDC) interconnections which provide increased controllability. In this con¬ 
text, we formulate a stochastic market-clearing algorithm that attains full spatio-temporal integration of 
reserve capacity, day-ahead and balancing markets. Against this benchmark we compare two deter¬ 
ministic market designs with varying degrees of coordination between the reserve capacity and energy 
services, both followed by a real-time mechanism. Our study reveals the inefficiency of determinis¬ 
tic approaches as the shares of wind power increase. Nevertheless, enforcing a tighter coordination 
between the reserves and energy trading floors may improve considerably the expected system cost 
compared to a sequential market design. Aiming to provide some insights for improvement of the se¬ 
quential market-clearing, we analyse the effect of explicit transmission allocation between energy and 
reserves for different HVDC capacities and identify the market dynamics that dictate the optimal ratio. 
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1 INTRODUCTION 


The European power system is undergoing significant restructuring both in terms of generation capa¬ 
bilities and market architecture. The high shares of variable and partially predictable generation, e.g., 
wind power, pose new challenges to system operators due to the increased reserve requirements aris¬ 
ing from the associated forecast errors. Now, although the appeal for improved efficiency and liquidity 
has already motivated the establishment of a single European day-ahead electricity market, the provi¬ 
sion and the activation of reserves remain mainly on a national level. However, in this new operation 
paradigm, the integration of the reserve capacity and balancing markets seems to be a necessary step 
in order to enable the cross-border exchange of regulating services and reduce the balancing costs. 

A major prerequisite for the deployment of cross-border reserves during real-time operation is the 
availability of transmission capacity between different areas of the power system. Meanwhile, the trans¬ 
mission grid should have adequate operational flexibility in order to allow for frequent re-dispatching 
of control reserves depending on the actual imbalances. HVDC interconnections can significantly 
contribute to a reliable solution owing to their increased controllability features that enable dynamic 
changes of the power flow direction which can be used to mitigate AC cross-border congestions Q]. 

Given the current European market setup, where reserve capacity, forward and balancing market- 
clearings are performed sequentially, the ex-ante allocation of transmission capacity between these 
trading floors is necessary. However, this procedure is subject to uncertainty due to the significant 
time interval between day-ahead scheduling and actual operation. The two main approaches to cope 
with this uncertainty are stochastic programming and deterministic models [2]. Despite the theoretical 
advances of the former, current operational practices are based on deterministic and usually static 
reserve requirements. For instance, in the new HVDC interconnection between Denmark and Norway 
(Skagerrak 4), 15% of its capacity will be permanently reserved for exchange of ancillary services [3]. 

Previous studies @], 0 have shown the benefit of cross-border reserve procurement in the North¬ 
ern European power system based on deterministic reserve requirements. In this work, we model 
and assess the efficiency of transmission capacity allocation for the exchange of balancing power in 
systems with a large stochastic wind power component. In this context, we define two deterministic 
market setups which are compared with the ideal stochastic clearing. In addition, we analyse the sen¬ 
sitivity of the deterministic solution with respect to particular parameters, i.e., wind power penetration 
and shares of transmission capacity between different markets. The proposed methodology is applied 
on a modified power system based on the IEEE RTS-96 with additional HVDC interconnections. 

2 MATHEMATICAL FORMULATION 

This section provides the mathematical formulation of three alternative market-clearing models em¬ 
ployed to settle the reserve capacity, day-ahead and balancing markets. All the market setups con¬ 
sidered in the present work are based on the following assumptions: /) only wind power uncertainty is 
considered assuming that load and availability of conventional generators are perfectly known, /'/) de¬ 
mand is completely inelastic and Hi) the market settles on an hourly basis and each trading period 
is independent since no inter-temporal constraints, e.g., ramping limits, are considered. The market¬ 
clearing models presented in the remaining of this section are described in a compact form, while the 
full model formulation is available in the Appendix. 

2.1 Stochastic energy and reserve co-optimization 

In order to determine the optimal amount of reserves and the day-ahead dispatch taking into account 
the anticipated real-time operation of the power system, market operator solves model|T]having a finite 
set of scenarios that captures the spatial dependence structure of the forecast errors. This formulation 
performs an intrinsic trade-off between the value of the available resources in the day-ahead (energy 
and reserve capacity) and the balancing (reserve deployment) stages, while it enables an implicit 
allocation of the transmission capacity among energy and reserve services. 
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where E s = {pc,Pw,r,8, fu} is the set of optimization variables and E^[-] the expec¬ 

tation operator over the scenario set which describes wind power uncertainty. This is a two-stage 
stochastic programming problem, where vectors pc and pw denote the day-ahead dispatch of con¬ 
ventional and wind power generators, respectively. Vector f contains the reserve capacity scheduled 
for units that provide balancing services. All re-dispatching actions, i.e., up- and down-regulation, wind 
spillage and load shedding, that take place in real-time operation are included in vector f w . The nodal 
voltage angles corresponding to day-ahead and balancing stage are denoted as 5 and 6 U , while W w 
corresponds to the realization of the random variable W, modelling the actual wind power produc¬ 
tion in scenario w. The power flows on the HVDC links, denoted as / at the day-ahead stage and 
f w during the balancing stage, are considered fully controllable and any voltage limitations and active 
power losses are neglected. Constants l, f and p w denote the load at each node of the system, the 
transmission capacity limits and the wind power installed capacity, respectively. 

Constraints (JTb) and (Qg} are the nodal power balance equations in the day-ahead and balancing 
stages. The day-ahead energy dispatch and reserve provision are limited by the generator maximum 
production through constraints {Tc} and (Qd). Similarly, constraints 07) ensure that the reserved ca¬ 
pacity of each generator does not exceed its reserve quantity offer R , while constraints dTFTb guarantee 
that the balancing actions are limited by the corresponding reserved capacity or rational constraints, 
i.e., load shedding and wind spillage are lower than or equal to the load and the realized wind power 
production, respectively. Finally, the set of constraints |TeJ and (HJ enforces the transmission capacity 
limits at the day-ahead and balancing stages. 

2.2 Deterministic market-clearing 

In order to avoid the computational burden of the stochastic market-clearing approach, market oper¬ 
ators employ deterministic models which have to procure reserve capacity based on explicit reserve 
requirements in order to deal with unforeseen events during real-time operation. This section provides 
the mathematical formulation of two deterministic approaches; the first allows for simultaneous pro¬ 
curement of energy and reserves while the second decouples completely these two services through 
a sequential settlement of the energy and reserve capacity markets. Both market organizations include 
a real-time balancing mechanism to compensate any deviations from the day-ahead schedule. 

2.2.1 Energy and reserve capacity co-optimization 

The day-ahead market-clearing problem for energy and reserve co-optimization is described through 
model [2l followed by the balancing market settlement of model [3] This market setup finds the opti¬ 
mal day-ahead dispatch (denoted by the superscript "*") and yields an implicit optimal transmission 
capacity allocation between energy and reserves, based on a single-valued forecast >V (conditional 
expectation) of wind power production, taking into account the explicit system-wide reserve require¬ 
ments RR S . The optimization variables of the day-ahead and the balancing market-clearing problems 
are included in the sets = {pc,Pw,r,S, /} and = {r,S,f}, respectively. 
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Day-ahead market-clearing 
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2.2.2 Sequential clearing of reserve capacity and energy markets 

The separation of reserve capacity and energy markets is implemented using models 0] and [5] This 
market architecture requires the explicit allocation of transmission capacity for reserve exchange, de¬ 
noted as X, as well as the definition of reserve requirements RR a for each area a of the power system. 
Model [^determines the least-cost reserve capacity f* which is used as input in the energy-only day- 
ahead clearing [5] whose optimization variables are contained in the set EfP = {pc,Pw,$, f}. Any 
corrective actions during real-time operation are compensated through a balancing market as of model 
[^considering that a single pool of control resources is accessible from all areas. 

Reserve capacity market-clearing Day-ahead market-clearing 
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3 RESERVE REQUIREMENTS AND WIND POWER UNCERTAINTY 

The reserve requirements are calculated using a probabilistic approach similar to 0, aiming to cover 
a pre-specified interval £ of the wind power predictive distribution. Write a,a € [0,1] the appropriate 
nominal proportions of the predictive cumulative distribution function (CDF) F such that £ = a —a, e.g., 
for reliability level £ = 99% it follows that a = 0.995 and a = 0.005. Then, the reserve requirements for 
up- and down-regulation from wind farm s, denoted as RRf and RR~ respectively, are calculated as 

RRf = (As - F- 1 (a)) f Wl RRJ = - p s ) fw (6) 

where ji s is the mean of F s . System-wide reserve requirements are obtained from ©, replacing F s 
with the CDF of the whole wind power portfolio and p s w with the total installed wind power capacity. 

The stochastic output of wind farm s, normalized by its nominal capacity, is modelled using a 
Beta distribution B s (a s ,/3 s ) as proposed in Q. To obtain a set Q of spatially correlated scenarios 
we follow the methodology presented in 0. Denote as X* the realization drawn from a multivariate 
Normal distribution A7(0, £), where £ the covariance matrix with diagonal elements equal to 1 and 
off-diagonal elements equal to the spatial correlation coefficient p of the corresponding wind locations. 
Then, the TV^ wind power production scenario is generated using the following transformation: 

yj = S(x-) —► W* = F~l (YJ) (7) 

where and F B j s denote the Gaussian and Beta CDF, respectively. 


4 


4 CASE STUDY 


In this section we evaluate the performance of the market-clearing algorithms described in section 
[2] using a power system comprised of areas 1 and 2 of the IEEE RTS-96 system presented in [9]. 
The offer price of each conventional generator in the day-ahead and balancing markets is set equal 
to the average of the three marginal cost segments provided in [9], while wind power offers its energy 
production at zero price. In addition, the minimum power production of all generators is assumed to 
be zero and the value of lost load is set to 1000$/MWh. OCGT plants provide their full capacity in 
the reserve market at cost per MW equal to 25% of their day-ahead offer price. The IGCC and CCGT 
units offer 40% and 25% of their capacity for reserve provision at a price equal to 10% and 5% of their 
marginal cost. Considering their limited flexibility, nuclear and coal units do not provide any balancing 
services. In order to focus our analysis on the HVDC tie-line, we follow a zonal pricing scheme similar 
to the NordPool market, neglecting the intra-area network constraints. This allows to model wind power 
production in areas 1 and 2 by employing two Beta distributions with shape parameters (a, /3) equal to 
(3.78, 1.62) and (5.67, 6.48), respectively. A set of 100 equiprobable scenarios is obtained assuming 
a correlation coefficient p equal to 0.35. 

Figure [T] shows the evolution of the expected power system operation cost as wind penetration 
increases, for the three considered market setups and an HVDC tie-line capacity of 200 MW. The ratio 
of the installed wind power capacity between areas 1 and 2 is fixed to 2:1 and wind power penetration 
is defined as a percentage of total system load. It can be observed that stochastic market-clearing out¬ 
performs both deterministic approaches in the whole range of wind penetration levels, exploiting the 
advanced information on the spatial characteristics of wind forecast uncertainty and managing to fully 
capture the benefits of cost-free wind power. The performance of the two deterministic market setups 
is similar up to a penetration level of 45%, while beyond this wind power share the deterministic energy 
and reserve co-optimization model achieves significantly lower expected costs. However, the ineffi¬ 
ciency of both deterministic models becomes apparent for large penetration levels, where increasing 
installed wind power capacity leads to higher expected cost. It should be noted that the curves in figure 
[Hare calculated using in-sample analysis, i.e., the same scenario set O describing wind uncertainty in 
the day-ahead stage is also used to calculate the expected balancing cost of each market design. 

In order to study the effect of the explicit transmission allocation X on the expected system cost, 
we consider a more penalizing reserve market where OCGT, IGCC and CCGT units provide 50%, 
40% and 40% of their capacity for reserve provision at a price equal to 15%, 30% and 30% of their 
day-ahead offer, respectively. For a constant wind penetration level of 24%, we perform a ’grid-search’ 
on the HVDC tie-line capacity and the parameter X, to obtain the corresponding expected system 
cost displayed in figure [2] The shape of this surface allows to identify three main directions aiming 
to analyse the effect of these parameters on the overall market efficiency. Note that the dotted line 
C represents the locus of the least-cost points for different HVDC capacities. Moving along direction 
A on the left-hand side of C, the expected cost reduces as the tie-line capacity increases since a 
larger pool of common resources is accessible from both areas. A similar trend applies for larger 
values of X (direction B ), where the reserve capacity market does not alter significantly the optimal 
day-ahead settlement by imposing additional constraints. For example, in a segmented market of 
reserve resources (low X), expensive units may be assigned to provide downward reserve, enforcing 
a corresponding lower bound on the energy dispatch, out of the merit order. Large values of X on 
high HVDC capacity (direction C ), allow the reserve market to pick more economical resources based 
solely on their capacity payment. However, this translates into increased expected balancing cost if 
mainly excess production situations occur during actual operation, i.e., generators with low marginal 
(and reserve capacity) cost are willing to buy back their production in lower price. 

5 CONCLUSIONS AND FUTURE WORK 

This paper considers the problem of optimal allocation of flexible HVDC transmission for the cross- 
border exchange of balancing reserves in presence of uncertain wind power generation. Setting as 
benchmark the stochastic energy and reserves co-optimization, we analyse the efficiency of two deter¬ 
ministic market designs to accommodate large shares of wind power and we show that both versions 
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Figure 1: Expected cost as a function of wind power penetration 
for different market-clearing models (X = 15%). 


Figure 2: Impact of FIVDC capacity and transmission allocation X on 
the expected cost of sequential market-clearing. 


behave fairly similar, only truly revealing the benefits from co-optimization when getting close to the 
system limits. Despite that the relative difference in market efficiencies may be highly sensitive to the 
system properties, the ranking is expected to stay the same. A wide range of plausible system setups 
is explored in order to perceive the market dynamics defining the optimal transmission allocation in the 
sequential market-clearing. It can be concluded that there generally exists an optimal allocation to be 
made, which however may dynamically vary depending on generation, load and system uncertainties. 

Future work may focus on the reformulation of the deterministic energy and reserve co-optimization 
algorithm by defining per area reserve requirements in order to reveal the optimal transmission capacity 
allocation of this model. This could serve as an insight towards an explicit rule that determines the 
optimal percentage of the tie-line capacity to be designated to reserve exchange in the sequential 
design. Further extensions could investigate the issues related to HVDC and AC network coordination 
as well as integrate the aspect of cross-border reserve provision in the HVDC expansion planning. 

BIBLIOGRAPHY 

[1] S. Chatzivasileiadis, “Power system planning and operation methods integrating the controllability 
of HVDC,” Ph.D. dissertation, ETH Zurich, 2013. 

[2] J. M. Morales, A. J. Conejo, H. Madsen, P. Pinson, and M. Zugno, Integrating Renewables in 
Electricity Markets: Operational Problems. Springer Science & Business Media, 2014. 

[3] Energinet.dk. (2011) Energinet.dk’s ancillary services strategy. [Online]. Available: 
https://www.energinet.dk/ 

[4] Y. Gebrekiros and G. Doorman, “Optimal transmission capacity allocation for cross-border ex¬ 
change of frequency restoration reserves (FRR),” in 18th Power Syst. Comput. Conf., 2014. 

[5] H. Farahmand and G. Doorman, “Balancing market integration in the northern european continent,” 
Applied Energy, vol. 96, pp. 316-326, 2012. 

[6] R. Doherty and M. O’Malley, “A new approach to quantify reserve demand in systems with signifi¬ 
cant installed wind capacity,” IEEE Trans. Power Syst., vol. 20, no. 2, pp. 587-595, 2005. 

[7] A. Fabbri, T. G. San Roman, J. R. Abbad, and V. H. M. Mendez Quezada, “Assessment of the cost 
associated with wind generation prediction errors in a liberalized electricity market,” IEEE Trans. 
Power Syst., vol. 20, no. 3, pp. 1440-1446, 2005. 

[8] S. Delikaraoglou and P. Pinson, “High-quality wind power scenario forecasts for decision-making 
under uncertainty in power systems,” in 13th Wind Integration Workshop, 2014. 

[9] H. Pandzic, Y. Dvorkin, T. Qiu, Y. Wang, and D. Kirschen. (2014) Library of the renewable 
energy analysis lab (REAL), University of Washington, Seattle, USA. [Online]. Available: 
http://www.ee.washington.edu/research/real/library.html 


6 












APPENDIX 


Nomenclature 

Indices and Sets 

n Index of system buses. 

i Index of dispatchable units. 

j Index of wind power units. 

I Index of transmission lines. 

uj Index of wind power scenarios. 

L ac Set of AC transmission lines. 

L dc Set of HVDC transmission lines. 

Ad/ Set of wind power units located at bus n. 
Ad/ Set of dispatchable units located at bus n. 


Parameters 


An 


Line-to-bus incidence matrix. The elements of the incidence matrix A, 
are equal to -1/1 if line l leaves/enters bus n and 0 otherwise. 

D n 


Demand at node n [MW]. 

Ci 


Day-ahead offer price of unit i [$/MWh]. 


r>R— 

Up/down reserve capacity offer price of unit i [$/MW]. 

C sh 

pW,max 

j 

Value of involuntarily shed load [$/MWh], 

pmax 
r i i 

Capacity of units i and j [MW]. 

dW 

j 


Forecast mean of wind power production [MW]. 

p* 

r ju) 


Wind power realization in scenario uj [MW]. 

r>+,max t )—,max 

^i ’ ■K'i 

Maximum up/down reserve provided by unit i [MW]. 

RR S ’ + 

, RR S ~ 

Maximum up/down system reserve requirements [MW]. 

/’max 

Jl 


Maximum capacity of line l [MW]. 

Bi 


Absolute value of the susceptance of AC line l [per unit]. 


Variables 


Pi,pf 

Rf,Rj 

r~y r. 

iuj ’ ICO 

r/ pi11 

"iuj 

/Sh 

l nuj 

Ki 


fAC fDC 
Jl > Jl 
fAC pC 

J luj ? J 1(jJ 


Day-ahead dispatch of units i and j [MW]. 

Up/down reserve capacity scheduled for unit i [MW]. 

Up/down regulation provided by unit i in scenario uj [MW]. 

Wind power spilled by unit j in scenario cu,[MW]. 

Load shedding at node n in scenario uj [MW]. 

Voltage angle at node n at the day-ahead market stage [rad]. 
Voltage angle at node n in scenario uj [rad]. 

Power flow in AC/HVDC line l at day-ahead market stage [MW]. 
Power flow in AC/HVDC line l in scenario uj [MW]. 
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Stochastic energy and reserve co-optimization 
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Deterministic energy and reserve capacity co-optimization 

Day-ahead market-clearing: Energy and reserve capacity 
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Balancing market-clearing 
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The set of optimization variables = {r+, r i ,p® pi11 , i® h , 8 ni f AC , fi DC , Vi, j, n, l}. 
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Sequential clearing of reserve capacity and energy markets 


Additional nomenclature: 


a, b 

M L ab 

Xi 

RR+,RR- 


E>+ D — 

ia > ia 


Index of system areas. 

Set of units located in area a. 

Set of lines connecting area a with area b. 

Transmission capacity of line l allocated to reserve exchange [%]. 
Up/down reserve requirements in area a [MW]. 

Up/down reserve capacity provided by unit i 
to cover requirements of area a [MW]. 


Reserve capacity market-clearing 
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Figure 3: Illustrative example: Capacity allocation for reserve exchange for a = {ai}, b = {a 2 , a 3 } - constraints |11g) and 111 hi . 
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Day-ahead market-clearing: Energy only 


S Y, c iPi 

u iei 
subject to 

E p7 + E n - A, - E A ‘n (/, AC + f, DC ) = 0. Vn 

j&Mi i&M ! n leL AC UL DC 

Pl <pr x -J2 R ta, Vi 

a 

Y, R ia<Pii Vi 


pf<Pj\ Vj 

fl AC = BlJ2 A lJn, V/ a 


?w 


AC 


- (i - Xj)fr x < f] Ac < a - ^)/r x , v/ € l ac 

- (i - ^)/r x < fi DC < a - ^)/r x , 

<5, = 0 

K > 0, Vi; p^ > 0, Vj; 

<5 n free, Vn; ff AC free, VI € L AC ] ff 0 free, V/ € L DC 


The set of optimization variables Q DA = {pi,pY ,S n , f AC , f t DC , Vz,j, n,/}. 


(12a) 


(12b) 

(12c) 

(12d) 

(12e) 

(12f) 

(I2g) 

(12h) 

( 121 ) 

(12j) 
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Balancing market-clearing 


Mim ) 

+ E ° sh C 

(13a) 

u iei 

nGAT 


subject to 



E (V-O+i? 

+ E C;-j ’T-pT") 


i&M ! n 



+ E A ln (ft C -fr + fP C -f? C )= 0, Vn 

(13b) 

leL AC UL DC 



f AC = B t E A ln 8 n , 

VI € L AC 

(13c) 


n 


a 

(13d) 

r~<J2 R ^ 

(13e) 

- // max < fi AC < // max , VI € 

(13f) 

- / z max < /^ C < ./) max , V/ € L dc 

(13g) 

Vn 

(13h) 

/>f : < rj- Vj 

(13i) 

II 

o 

(13j) 

rf^i >0 , Vi; pf" > Q.Vj: C>0 ,Vn; 

<5 n free, Vn; f AC free, V/ € L AC ; /j DC free, V/ € 

(13k) 


The set of optimization variables ® B = {rf ,r i ■ ^ h , 8 n , f AC , , Vi, j, n, l}. 

Constraints (13c) and (13d) imply that reserve deployment depends only on reserve availability 
in real-time, irrespective of the initial area requirements. There is a single pool of balancing 
resources accessible from all areas during real-time operation. 
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